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ABSTRACT 

Earthquake  data  from  California  from  the  period  of 
January  1,  1969,  to  December  31,  1971,  has  been  analyzed 
by  two  methods.   The  first  method  compares  the  occurrence 
of  earthquakes  with  the  Synodic,  Draconic  and  Anomalistic 
Lunar  Periods.   Statistical  tests  to  determine  if  the  earth- 
quakes are  uniformly  distributed  over  each  of  the  three 
lunar  periods  are  presented.   The  second  method  compares 
three  components  of  the  tidal  force  vector  and  their 
derivatives  at  the  times  of  earthquake  occurrences. 


TABLE  OF  CONTENTS 

I.   INTRODUCTION  -------------------   5 

II.   ANALYSIS  AND  DISCUSSION  OF  THE  LUNAR  PERIOD 

COMPUTER  PROGRAM  -----------------  13 

III.   RESULTS-  - _  -  _  -  _  -  -  -  17 

A.  SAN  FERNANDO  VALLEY  EARTHQUAKE  DATA-  -----  19 

1.  Synodic  Month-  ---------  19 

2.  Draconic  Month  --------------19 

3.  Anomalistic  Month-  ------------20 

B.  SAN  FRANCISCO  BAY  EARTHQUAKE  DATA  I-  -----  24 

1.  Synodic  Month-  --------------    2^ 

2.  Draconic  Month  ---- ---26 

3.  Anomalistic  Month-  ------------28 

C.  SAN  FRANCISCO  BAY  EARTHQUAKE  DATA  II  -----  30 

1.  Synodic  Month-  --------------  30 

2.  Draconic  Month  --------------32 

3.  Anomalistic  Month-  ------------  34 

D.  TIDAL  FORCE  COMPONENTS  WITH  RESPECT  TO 

SAN  FERNANDO  VALLEY  DATA  -----------  36 

IV.   DISCUSSION  OF  RESULTS  AND  CONCLUSIONS-  ------  42 

V.   THE  COMPUTER  PROGRAM  ---------------  45 

LIST  OF  REFERENCES  ----- _  _  -  -  -  51 

INITIAL  DISTRIBUTION  LIST-  ----------  53 

FORM  DD  1473  - __-__-- 54 


ACKNOWLEDGEMENTS 

I  wish  to  express  my  deep  appreciation  to  Professor  Rex 
Shudde  and  Professor  Donald  Barr,  whose  valuable  criticism 
and  suggestions  contributed  substantially  to  the  completion 
of  this  thesis.   I  would  also  like  to  thank  M.  D.  Wood,  Ph.D., 
of  the  National  Center  for  Earthquake  Research,  Menlo  Park, 
California,  for  providing  the  data  used  in  this  thesis  as 
well  as  the  use  of  his  extensive  reference  literature  per- 
taining to  earthquakes.   Further,  I  would  like  to  extend 
my  appreciation  to  the  staff  at  the  W.  R.  Church  Computer 
Center,  Naval  Postgraduate  School  for  their  assistance. 


I.   INTRODUCTION 

A  lunar  correlation  with  the  time  of  the  occurrence 
of  earthquakes  as  reported  by  Allen  [Ref.  1]  has  been  one 
of  the  methods  used  in  trying  to  determine  the  causes  of 
earthquakes  for  the  better  part  of  a  century.   Knopoff 
[Ref.  3]  reported  that  it  may  be  that  earthquakes  have  a 
greater  tendency  to  be  triggered  when  the  sun  and  the  moon 
are  nearly  aligned  with  the  earth.   An  interesting  observa- 
tion which  was  reported  by  Wigand  [Ref.  16]  was  that  the 
sun  and  the  moon  were  almost  aligned  with  the  earth  at 
the  time  of  the  San  Fernando  Valley  Earthquake  of 
February  9,  1971.   It  is  further  contended  by  Wigand 
[Ref.  16]  that  there  is  strong  evidence  that  some  earth- 
quakes are  related  to  tidal  force,  and  some  are  actually 
triggered  by  the  tidal  force.   The  objective  of  this 
thesis  was  to  investigate  these  findings  by  analyzing  the 
gravitational  forces  of  the  moon  and  the  sun  with  relation 
to  recorded  earthquakes. 

Two  approaches  in  analyzing  earthquake  data  are 
presented.   The  first  approach  was  to  compare  the  occurrence 
of  earthquakes  with  the  Synodic  Lunar  Period,  the  Draconic 
Lunar  Period,  and  the  Anomalistic  Lunar  Period.   The  second 
approach  was  to  investigate  the  tidal  force  exerted  by  the 
sun  and  the  moon  with  relation  to  155  recorded  earthquakes 
which  occurred  in  the  San  Fernando  Valley  from  April,  1971, 
to  November,  1971. 


The  three  lunar  periods  discussed  in  the  first  approach 
are  as  follows: 

1.  Synodic  Lunar  Period,  mean  value  29.530589  solar 
days.   Defined  as  new  moon  to  new  moon,  where  the 
new  moon  occurs  when  the  geocentric  longitudes  of 
the  sun  and  the  moon  are  the  same  [Ref.  12]. 

2.  Draconic  Lunar  Period,  mean  value  27.212220  solar 
days .   Defined  as  the  ascending  node  to  the 
ascending  node,  where  the  ascending  node  pas- 
sage of  the  moon  is  from  minus  degrees  latitude 
to  positive  degrees  latitude  [Ref.  12]. 

3.  Anomalistic  Lunar  Period,  mean  value  27.554531 
solar  days.   Defined  as  perigee  to  perigee, 
where  perigee  is  that  point  of  the  moon's  orbit 

at  which  the  moon  is  nearest  to  the  earth  [Ref.  12] 
Five  lunar  periods  [Ref.  12]  actually  exist  but  the 
above  three  were  included  for  the  reasons  given  below. 
In  analyzing  the  earthquake  data  with  respect  to  the 
Synodic  Lunar  Period,  the  occurrence  of  earthquakes  with 
respect  to  the  positions  of  the  sun  and  the  moon  are 
examined.   In  analyzing  the  data  with  respect  to  the 
Draconic  Lunar  Period,  the  occurrence  of  earthquakes  with 
respect  to  the  moon's  position  in  the  Northern  Hemisphere 
and  Southern  Hemisphere  were  examined.   This  may  have  been 
important  since  the  data  used  included  only  earthquakes 
which  occurred  in  California.   In  analyzing  the  earthquake 
data  with  respect  to  the  Anomalistic  Lunar  Period,  the 


difference  in  the  gravitational  force  of  the  moon  on  the 
earth  was  examined  as  the  moon  proceeded  from  perigee  to 
apogee  to  perigee. 

The  earthquake  data  received  from  the  National  Center 
for  Earthquake  Research  spanned  three  time  periods  and 
occurred  in  two  areas  of  California.   The  first  set  of  data 
included  155  earthquakes  recorded  from  April  1,  1971,  to 
November  17,  1971.   These  are  classified  as  aftershocks  of 
the  February  9,  1971,  San  Fernando  Valley  Earthquake 
(Figure  1).   The  magnitudes  ranged  from  1.0  to  4.2  as 
measured  on  the  Richter  Scale.   This  set  of  data  will  be 
referred  to  as  the  San  Fernando  Valley  Data  (SFVD) . 

The  second  set  of  data  included  2879  earthquakes 
recorded  from  January  1,  1969,  to  August  23,  1970.   These 
earthquakes  occurred  in  the  San  Francisco  Bay  Area  (Figure 
2)  in  the  vicinity  of  the  San  Andreas  Fault,  the  Hayward 
Fault,  and  the  Calaveras  Fault.   The  magnitudes  ranged 
from  0.11  to  5.06  as  measured  on  the  Richter  Scale  [Refs. 
4  and  5].   This  set  of  data  will  be  referred  to  as  the  San 
Francisco  Bay  Area  Data  I  (SFBAD-I). 

The  third  set  of  data  included  2734  earthquakes  recorded 
from  August  1,  1970,  to  December  31,  1971.   The  area  was 
the  same  as  SFBAD-I.   The  magnitudes  ranged  from  0.29  to 
4.45  as  measured  on  the  Richter  Scale  [Refs.  5  and  6]. 
This  set  of  data  will  be  referred  to  as  the  San  Francisco 
Bay  Area  Data  II  (SFBAD-II). 
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FIGURE  2 


SFBAD-I  and  SFBAD-II  both  included  earthquakes  which 
occurred  in  August,  1970.  In  the  analysis,  as  discussed 
in  the  next  section,  the  duplicate  earthquakes  were  eliminated 

In  analyzing  the  data  with  respect  to  the  lunar  periods, 
the  following  hypotheses  were  investigated: 

Null  Hypothesis:   Earthquakes  occur  in  a 

uniform  distribution  with 
respect  to  the  lunar  period. 
Alternative  Hypothesis :   Earthquakes  do  not 

occur  in  a  uniform  distri- 
bution with  respect  to  the 
lunar  period. 
Two  tests  were  used  in  the  analysis:   The  Chi-Square 
Test  for  Significance  [Ref.  17],  and  the  Goodness  of  Fit 
Statistic  Vn  [Ref.  13] .   The  tests  are  discussed  in  what 
follows : 

The  Chi-Square  Test  for  Significance 

N,  the  total  number  of  observations  in  the 

samples  which  were  analyzed. 
H,  the  total  number  of  equal  intervals  in 
which  the  lunar  period  was  divided. 

X  =  T7  ,  the  expected  number  of  observations  \vhich 
occurred  in  each  interval  of  a  uniform 
distribution. 
Ti,  the  total  number  of  observations  recorded  within 
each  interval,  i=l,...,H. 
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df  =  H  -  1,  the  degrees  of  freedom  associated  with 
the  test. 

2      H   (T  -X)2 

X     .     =      T.      — =! ,  the  Chi-Square  statistic. 

df   j  =  i    X 

The  Goodness  of  Fit  Statistic  Vn : 

N,  the  total  number  of  observations  in  the 

sample  analyzed. 
F(x),  the  uniform  distribution  on  the  lunar 

period  in  which  the  origin  was  the  origin 

of  the  lunar  period. 
F^(x),  the  distribution  of  the  observations  on 

the  lunar  period  with  the  same  origin  as  F(x). 
X,  the  time  of  each  earthquake  in  the  sample. 
Vn  =  Sup  (Fj^(x)-F(x))  -  Inf  (F^  (x)  -  F  (x) )  , 

-oo<x<oo  _oo<x<oo 

,  the 
Vn  statistic. 

Stephens  [Ref.  13]  stated  that  if  the  observations 
are  points  on  a  circle,  the  value  of  Vn  obtained  from  the 
above  equation  does  not  depend  on  the  choice  of  the  origin 
for  measuring  x.   The  distribution  of  Vn  refers  to  the  dis- 
tribution under  the  null  hypothesis  [Ref.  13] .   Tables  of 
significance  of  the  statistic  /N  Vn  are  included  in 
Reference  13. 

A  simulation  of  1000  earthquakes  over  a  period  of  six 
months  was  included  in  the  analysis.   The  simulation  was 
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performed  to  test  the  uniform  distribution  of  the  events 
over  the  three  lunar  periods. 

The  final  analysis  in  this  thesis  v\[as  made  possible 
by  the  use  of  a  computer  program  by  Professor  Rex  H. 
Shudde  of  the  Naval  Postgraduate  School.   This  computer 
program  calculated  three  components  of  the  tidal  force 
(measured  in  microgals)  at  the  epicenter  of  an  earthquake. 

The  three  components  are  the  radial  component,  the 
north  component,  and  the  east  component.   The  radial  com- 
ponent of  the  tidal  force  vector  is  through  the  center  of 
the  earth  and  the  epicenter  of  the  earthquake.   The  north 
and  east  components  of  the  tidal  force  vector  are  in  the 
plane  which  passes  through  the  epicenter  and  is  perpendicular 
to  the  radial  component.   The  computer  program  also  calcu- 
lated the  rates  of  change  of  the  three  components  described 
above.   The  three  components  were  plotted  for  a  24-hour 
time  period  on  the  day  each  earthquake  occurred.   A  visual 
comparison  of  the  plots  of  the  155  earthquakes  which 
occurred  in  the  San  Fernando  Valley  with  respect  to  the 
three  components  of  the  tidal  force  was  made.   The  results 
of  the  comparison  are  in  Section  III. 
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II.   ANALYSIS  AND  DISCUSSION  OF  THE  LUNAR 
PERIOD  COMPUTER  PROGRAM 

The  computer  program  used  to  analyze  the  frequency  of 
earthquake  occurrences  during  the  Synodic,  Draconic  and 
Anomalistic  Lunar  Periods  contains  the  following  steps. 
The  first  step  involves  the  time  frame  of  the  three  sets 
of  data  and  the  lengths  of  the  lunar  periods  within  these 
time  frames.   The  mean  times  of  the  three  lunar  periods, 
as  given  in  the  Introduction,  were  not  used  in  the  program, 
since  the  actual  lengths  of  the  periods  vary  considerably. 
An  example  of  this  difference  was  the  Anomalistic  Month 
which  has  a  mean  period  of  27.32166  days.   During  1971  the 
length  of  the  Anomalistic  Month  ranged  from  25.41667 
days  to  28.45833  days  [Ref.  9]. 

The  times  of  new  moon  and  lunar  perigee  for  the  years 
1969,  1970  and  1971  are  tabulated  in  References  7,  8  and 
9.   From  these  times  (to  be  referred  to  as  origins)  the 
different  period  lengths  of  the  Synodic  and  Anomalistic 

onths  were  computed.   References  7,  8  and  9  also  include 
tables  indicating  the  apparent  latitude  of  the  moon  for 
the  respective  years  at  intervals  of  one  half  day.   The 
origins  for  the  Draconic  Month,  as  the  moon  passed  from  the 
Southern  Hemisphere  to  the  Northern  Hemisphere,  were 
computed  from  these  tables  by  linear  interpolation.   This 
procedure  should  be  fairly  accurate  since  it  involves  only 
minutes  and  seconds  of  latitude. 
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The  origins  of  the  three  lunar  periods  and  the  times  of 
each  earthquake  occurrence  were  put  into  the  program  in 
year,  month,  day,  hours,  minutes  and  seconds  in  Greenwich 
Mean  Time  (GMT).   Within  the  program  the  year,  month  and 
day  of  each  event  were  converted  to  the  respective  Julian 
Day  Number  (JD)  [Ref.  14].   The  hours,  minutes  and  seconds 
constituted  a  fraction  (Fract)  of  a  day.   The  Julian  Date 
"was  computed  as  follows: 

Julian  Date  =  JD  +  Fract  -  0.50 

The  one  half  day  was  subtracted  due  to  the  differences 
in  origins  of  the  Julian  Day  Number  and  GMT.   The  Julian 
Day  Number  had  its  origin  at  noon  whereas  GMT  has  its  origin 
at  midnight. 

In  order  to  simplify  the  statistical  tests  which  followed 
in  the  program,  the  data  was  truncated  to  insure  that  no 
fractions  of  a  period  were  included.   An  example  of  this 
was  the  SFBAD- I  and  the  Synodic  Month.   The  time  the 
SFBAD- I  started  was  January  1,  1969,  which  was  between  the 
two  new  moons  of  December  19,  1968,  and  January  18,  1969. 
Since  the  events  that  may  have  occurred  from  December  19, 
1968,  to  December  31,  1968,  were  not  included,  the  events 
from  January  1,  1969,  to  January  18,  1969,  were  excluded 
in  the  analysis.   Similarly  fractions  of  periods  were 
eliminated  for  the  Draconic  and  Anomalistic  Months. 

Each  of  the  three  lunar  time  periods  was  reduced  to  the 
unit  circle  as  explained  in  the  following  example.   The 
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length  of  the  Synodic  Lunar  Period  from  April  25,  1971, 
to  May  24,  1971,  was  29.354170  days.   At  the  time  of  the  new 
moon  on  April  25,  the  value  on  the  unit  circle  was  zero. 
At  the  time  of  the  first  quarter  moon  on  May  2,  the  value 
on  the  unit  circle  was  0.25.   At  the  time  of  the  full  moon 
on  May  9,  the  value  on  the  unit  circle  was  0.50.   At  the 
time  of  the  third  quarter  moon  on  May  17,  the  value  on  the 
unit  circle  was  0.75.   Finally  at  the  time  of  the  new  moon 
on  May  24,  the  value  on  the  unit  circle  was  1.0  [Ref.  9]. 

The  time  at  which  each  earthquake  occurred  during  the 
lunar  period  (residual)  was  placed  upon  the  unit  circle 
with  respect  to  the  respective  origins.   This  was  accom- 
plished by  dividing  the  difference  between  the  Julian  Date 
of  the  earthquake  and  the  Julian  Date  of  the  origin  immedi- 
ately preceding  the  earthquake  by  the  length  of  the  period 
in  which  the  earthquake  occurred. 

The  unit  circle  was  divided  into  20  equal  intervals 
in  order  to  perform  the  Chi-Square  Test  for  Significance 
[Ref.  17]  as  explained  in  the  Introduction.   Ostle  [Ref.  10] 
suggests,  as  a  rule  of  thumb,  to  use  no  fewer  than  three  for 
the  expected  frequency  in  each  interval.   Ostle  [Ref.  10] 
further  contends  that  if  some  expected  numbers  are  too  small 
the  Chi-Square  statistic  will  be  a  poor  indicator  of  the 
validity  of  the  hypothesis  under  test.   By  selecting  20 
intervals,  the  earthquakes  in  the  three  sets  of  data  comply 
with  the  rule  of  thumb  for  expected  frequencies. 
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The  residuals  that  were  computed  were  then  sorted  into 
ascending  order  [Ref.  2]  so  that  the  Vn  Goodness  of  Fit 
Test  could  be  computed  as  explained  in  the  Introduction. 
The  Vn  Goodness  of  Fit  Test  is  a  modification  of  the 
Kolmogorov- Smirnov  Test  [Ref.  13]  which  is  considered  to 
be  generally  a  more  powerful  test  than  the  Chi-Square 
Test  [Ref.  10]. 

The  output  of  the  program  consisted  of  the  total  number 
of  events  that  occurred  during  each  of  the  three  lunar 
periods;  the  number  of  events  and  the  percentages  whicli 
occurred  within  each  of  the  20  intervals,  the  Chi-Square 
statistic,  and  the  Vn  and  /M  Vn  statistics. 
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III.   RESULTS 

The  simulation  of  the  1000  earthquakes  as  stated  in  the 
Introduction  was  performed  to  test  the  uniform  distribution 
of  the  events  on  the  three  lunar  periods.   The  null  hypothe- 
sis that  the  simulated  earthquakes  were  distributed  uniform- 
ly was  accepted  for  each  of  the  lunar  periods.   The  bases 
for  the  acceptances  were  the  Chi-Square  Test  with  19 
degrees  of  freedom  and  the  Vn  Test.   The  Chi-Square  statis- 
tics and  the  /N  Vn  statistics  calculated  are  given  below: 

Synodic   Draconic   Anomalistic 
.  Chi-Square      23.28      22.04       28.96 

/N  Vn  1.4002     1.2130      1.3995 

The  null  hypothesis  and  alternative  hypothesis  are: 
Null  Hypothesis:   Earthquakes  occur  in  a 

uniform  distribution  with  respect 
to  the  lunar  period. 
Alternative  Flypothesis:   Earthquakes  do  not  occur 

in  a  uniform  distribution  with 
respect  to  the  lunar  period. 
If  the  null  hypothesis  is  accepted,  it  might  imply  that 
the  moon's  gravitational  effects  on  the  earth  do  not  trigger 
earthquakes.   Conversely  if  the  null  hypothesis  is  rejected 
and  the  alternative  hypothesis  is  accepted  then  it  might 
imply  that  the  moon's  gravitational  effects  on  the  earth 
do  trigger  earthquakes. 
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The  results  of  the  statistical  tests  performed  on  the 
three  sets  of  data  with  respect  to  the  Synodic,  Draconic 
and  Anomalistic  Lunar  Periods  are  presented  in  tables. 
Results  suggested  by  the  tidal  force  plots  are  discussed 
following  the  tables. 

In  each  table  the  following  definitions  are  used: 
Mag  -  Magnitudes  of  earthquakes  greater  than  or 

equal  to  the  magnitude  entered  in  the  table. 
No   -  Number  of  earthquakes  in  the  sample. 
Chi  -  The  Chi-Square  statistic  with  19  degrees 

of  freedom. 
a    -  The  one  per  cent  significance  level  for  the 

Chi-Square  statistic  [Ref.  11]. 
Vn   -  The  Vn  statistic. 

Vn/N  -  The  Vn  statistic  for  large  values  of  No. 
b    -  The  one  per  cent  significance  level  for  the  Vn 

statistic . 
Ace  -  Acceptance  of  the  null  hypothesis  at  the  one 

per  cent  significance  level. 
Rej  -  Rejection  of  the  null  hypothesis  at  the  one 

per  cent  significance  level. 
Polar  plots  of  the  occurrence  of  earthquakes  during  the 
lunar  periods  are  presented  in  Figures  3,  4,  5,  6,  7,  8,  9, 
10,  and  11. 
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A.   SAN  FERNANDO  VALLEY 

1,  Synodic  Month 

Mag    No    Chi  a  Vn  Vnv^      b 

All    110    43.818  REJ  .18258  1.91496  ACC 

>  1.75   106    33.245  ACC  .15788  1.62552  ACC 

>  2.00   104    34.077  ACC  .15444  1.57496  ACC 

>  2.25    96    38.167  REJ  .14963  1.46607  ACC 

>  2.50    88    28.364  ACC  .14395  1.35035  ACC 

Table  1 

2.  Draconic  Month 

Mag    No    Chi  a  Vn  Vn/N      b 

>  All    119    34.109  ACC  .16705  1.82225  ACC 

>  1.75   115    26.739  ACC  .14775  1.58449  ACC 

>  2.00   112    25.857  ACC  .14285  1.51158  ACC 

>  2.25   103    28.068  ACC  .14267  1.44793  ACC 

>  2.50    93    26.355  ACC  .14475  1.39593  ACC 

Table  2 
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3.   Anomalistic  Month 

Mag  No    Chi     a      Vn     Vn/N     b 

All  122  21.934  ACC  .15824  1.74781  ACC 

>  1.75  118  19.627  ACC  .14189  1.54135  ACC 

>  2.00  115  20.478  ACC  .13843  1.48455  ACC 
>_   2.25  104  21.000  ACC  .14804  1.50976  ACC 

>  2.50  93  19.903  ACC  .14287  1.37782  ACC 


Table  3 
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B.   SAN  FRANCISCO  BAY  EARTHQUAKE  DATA  I 

1.   Synodic  Month 

Mag  No  Chi  a     Vn  Vn/H  b 

All  2721  156.405  REJ  .05351  2.79130  REJ 

>  0.50  2554  171.090  REJ  .05848  2.95563  REJ 

>  0.75  2333  187.660  REJ  .06551  3.16441  REJ 

>  1.00  2008  202.358  REJ  .06693  2.99933  REJ 

>  1.25  1580  204.278  REJ  .07694  3.03841  REJ 

>  1.50  1193  186.966  REJ  .08954  3.09260  REJ 

>  1.75  803  170.425  REJ  .09327  2.64308  REJ 

>  2.00  483  106.772  REJ  .09356  2.05609  REJ 

>  2.25  301  93.485  REJ  .11281  1.95726  ACC 

>  2.50  177  125.034  REJ  .16622  2.21147  REJ 

>  2.75  111  98.189  REJ  .18296  1.92756  ACC 

>  3.00  64  83.500  REJ  .24541  1.96325  REJ 

Table  4 
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2. 

Draconic  Month 

Mag 

No 

Chi 

a 

Vn 

Vn/N 

b 

All 

2780 

127.611 

REJ 

.04939 

2, 

.60438 

REJ 

> 

0.50 

2611 

141.823 

REJ 

.04908 

2. 

.50802 

REJ 

> 

0.75 

2384 

162.342 

REJ 

.05230 

2, 

.55383 

REJ 

> 

1.00 

2048 

163.230 

REJ 

.06090 

2, 

.75592 

REJ 

> 

1.25 

1611 

180.521 

REJ 

.06805 

2, 

.73115 

REJ 

> 

1.50 

1214 

148.965 

REJ 

.07410 

2, 

.58155 

REJ 

> 

1.75 

818 

115.203 

REJ 

.08995 

2, 

.57271 

REJ 

> 

2.00 

493 

65.823 

REJ 

.09078 

2, 

.01568 

REJ 

> 

2.25 

310 

48.193 

REJ 

.11873 

2, 

.09050 

REJ 

> 

2.50 

184 

70.348 

REJ 

.14956 

2, 

.02866 

REJ 

> 

2.75 

116 

58.138 

REJ 

.13858 

1, 

.49258 

ACC 

> 

3.00 

66 

45.515 
Table  5 

REJ 

.14618 

1, 

.18757 

ACC 
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3.   Anomalistic  Month 

Mag  No      Chi     a     Vn  Vn/N    b 

All  2794  148.190  REJ  .05336  2.82060  REJ 

>  0.50  2625  144.379  REJ  .05320  2.72548  REJ 
_>  0.75  2398  153.409  REJ  .05437  2.66237  REJ 

>  1.00  2061  156.613  REJ  .06303  2.86128  REJ 

>  1.25  1623  142.878  REJ  .06268  2.52532  REJ 

>  1.50  1224  115.314  REJ  .07092  2.48117  REJ 

>  1.75  822     94.448  REJ  .08887  2.54783  REJ 

>  2.00  496     51.823  REJ  .09322  2.07606  REJ 

>  2.25  310     40.710  REJ  .10603  1.86686  ACC 

>  2.50  184     38.826  REJ  .12593  1.70822  ACC 

>  2.75  116     34.345  ACC  .17119  1.84375  ACC 

>  3.00  66     26.727  ACC  .14634  1.18887  ACC 

Table  6 
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C.   SAN  FRANCISCO  BAY  EARTHQUAKE  DATA  II 
1.   Synodic  Month 

Mag      No  Chi     a     Vn  Vn/N    b 

All     2595  135.767  REJ  .06405  3.26256  REJ 

>_   0.50     2571  132.034  REJ  .06204  3.14580  REJ 

>  0.75     2438  132.336  REJ  .06392  3.15630  REJ 

>  1.00     2144  122.847  REJ  .06422  2.97360  REJ 

>  1.25     1691  96.002  REJ  .06416  2.63846  REJ 

>  1.50     1228  78.384  REJ  .07451  2.61101  REJ 

>  1.75      810  65.506  REJ  .08914  2.53705  REJ 

>  2.00      459  47.710  REJ  .11265  2.41349  REJ 
>_  2.25     275  35.764  ACC  .11716  1.94280  ACC 

>  2.50      155  33.774  ACC  .14005  1.74338  ACC 
_>  2.75      79  27.076  ACC  .18849  1.67538  ACC 

Table  7 
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2,   Draconic  Month 


Mag 

No 

Chi 

All 

2556 

31.512 

> 

0.50 

2534 

32.140 

> 

0.75 

2406 

31.506 

> 

1.00 

2117 

23.680 

> 

1.25 

1668 

17.875 

> 

1.50 

1212 

25.294 

> 

1.75 

810 

18.151 

> 

2.00 

454 

26.529 

> 

2.25 

271 

26.491 

> 

2.50 

152 

31.158 

> 

2.75 

77 

25.597 
Table  8 

a     Vn  Vn/N    b 

ACC  .03189  1.61201  ACC 

ACC  .03155  1.58830  ACC 

ACC  .03599  1.76540  ACC 

ACC  .03257  1.49860  ACC 

ACC  .03807  1.55467  ACC 

ACC  .05430  1.89038  ACC 

ACC  .05296  1.49892  ACC 

ACC  .06228  1.32712  ACC 

ACC  .07411  1.22006  ACC 

ACC  .09405  1.15959  ACC 

ACC  .13345  1.17102  ACC 
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3.   Anomalistic  Month 

Mag  No      Chi     a     Vn  Vn/N  b 

All  2634  36.752  REJ  .04930  2.53041  REJ 

>  0.50  2611  35.181  ACC  .04856  2.48153  REJ 

>  0.75  2479  40.943  REJ  .05378  2.67761  REJ 

>  1.00  2184  45.982  REJ  .05309  2.48120  REJ 

>  1.25  1717  43.641  REJ  .05110  2.11722  REJ 

>  1.50  1251  48.265  REJ  .06704  2.37127  REJ 

>  1.75  829  49.673  REJ  .07730  2.22563  REJ 

>  2.00  473  36.302  REJ  .09228  2.00702  REJ 

>  2.25  286  31.203  ACC  .10886  1.84091  ACC 

>  2.50  162  32.321  ACC  .11470  1.45990  ACC 
>_  2.75  84  31.238  ACC  .14562  1.33464  ACC 

Table  9 
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D.   TIDAL  FORCE  COMPONENTS  WITH  RESPECT  TO  THE 
SAN  FERNANDO  VALLEY  DATA 

Figure  12  is  an  illustration  of  the  tidal  force  plots 
used  in  the  analysis.   The  occurrence  of  the  earthquake 
in  the  plot  was  May  19,  1971  at  0451  GMT.   The  magnitude 
of  the  earthquake  was  2.6.   The  radial  component  of  the 
tidal  force,  as  measured  in  microgals,  was  positive  and 
decreasing  in  value  (negative  slope).   The  north  component 
was  negative  and  increasing  in  value  (positive  slope) . 
The  east  component  was  negative  and  decreasing  in  value 
(negative  slope) . 

The  same  110  SFVD  earthquakes  which  were  analyzed  with 
respect  to  the  Synodic  Lunar  Period  were  investigated  as 
the  event  presented  above.   The  results  are  tabulated  in 
Table  10.   The  time  frame  in  which  these  earthquakes 
occurred  was  from  April  25,  1971,  to  October  19,  1971. 

Ten  thousand  earthquakes  uniformly  distributed  in  time 
were  simulated  during  the  previously  mentioned  time  frame 
and  in  the  same  area  as  the  SFVD.   These  simulated  earth- 
quakes were  analyzed  in  the  same  manner.   The  results  are 
also  tabulated  in  Table  10. 

Each  component  at  the  time  of  the  earthquake  was  either 
positive  or  negative  in  value.   The  derivatives  of  each 
of  the  three  components  was  either  positive  or  negative 
at  the  time  of  the  earthquake.   Sixty-four  combinations  of 
the  three  components  and  their  slopes  were  then  considered. 
In  Table  10  the  follovv'ing  definitions  are  used. 
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RAD   -   The  radial  component  of  the  tidal  force. 
NOR   -   The  north  component  of  the  tidal  force. 
EAS   -   The  east  component  of  the  tidal  force. 
REC   -   The  number  of  earthquakes  recorded  with  the 

same  combination. 
PREC  -   The  percentage  of  the  earthquakes  recorded 

with  the  same  combination. 
SIM   -   The  number  of  simulated  earthquakes  recorded 

with  the  same  combination. 
PSIM  -   The  percentage  of  the  simulated  earthquakes 

with  the  same  combination. 
+    -   Component  of  tidal  force  was  positive. 

Component  of  tidal  force  was  negative. 
p  -  Value  of  the  component  was  increasing. 
n    -   Value  of  the  component  was  decreasing. 

The  data  in  Table  10  were  used  to  compare  the  distribu- 
tion of  the  earthquakes  with  the  distribution  of  the  simulated 
earthquakes  by  means  of  a  Contingency  Table  [Ref .  10] . 
The  null  hypothesis  was  that  the  two  distributions  are  the 
same.   The  combinations  of  the  three  components  of  the 
tidal  force  vector  and  the  three  derivatives,  for  which 
the  actual  and  the  simulated  earthquakes  never  occurred, 
were  eliminated.   This  resulted  in  a  2X34  Contingency  Table. 
The  Chi-Square  statistic  with  53  degrees  of  freedom  was 
computed  (39.972)  and  the  null  hypothesis  was  accepted. 
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RAD      NOR      EAS       REC      PREC      SIM      PSIM 


+p 

+P 

+P 

0 

0 

0 

0 

-p 

+P 

+P 

5 

4.55 

206 

2.06 

+n 

+  P 

+P 

0 

0 

0 

0 

-n 

+P 

+P 

0 

0 

0 

0 

+P 

-P 

+P 

0 

0 

21 

.21 

-P 

-P 

+P 

4 

3.64 

469 

4.69 

+n 

-P 

+P 

0 

0 

0 

0 

-n 

-P 

+P 

0 

0 

1 

.01 

+P 

+  n 

+P 

0 

0 

74 

.74 

-P 

+n 

+P 

3 

2.73 

179 

1.79 

+n 

+  n 

+P 

0 

0 

0 

0 

-n 

+n 

+P 

0 

0 

0 

0 

+P 

-n 

■^P 

1 

.91 

215 

2.15 

-P 

-n 

+P 

9 

8.18 

1350 

13.50 

+n 

-n 

+P 

0 

0 

0 

0 

-n 

-n 

+P 

0 

0 

6 

.06 

+P 

+P 

-P 

0 

0 

0 

0 

-P 

+P 

-P 

0 

0 

0 

0 

+n 

+P 

-P 

1 

.91 

69 

.69 

-n 

+P 

-P 

5 

4.55 

161 

1.61 

+P 

-P 

-P 

0 

0 

0 

0 

-P 

-P 

-P 

0 

0 

17 

.17 

+n 

-P 

-P 

4 

3.63 

246 

2.46 

-n 

-P 

-P 

14 

12.73 

1356 

13.56 

+p 

+n 

-P 

0 

Table  10 

0 

0 

0 
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RAD 
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PREC 

SIM 
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-P 

+n 

-P 

0 

0 

0 

0 

+n 

+n 

-P 

0 

0 

0 

0 

-n 

+n 

-P 

1 

.91 

236 

2.36 

+P 

-n 

-P 

0 

0 

0 

0 

-P 

-n 

-P 

1 

.91 

5 

.05 

+n 

-n 

-P 

0 

0 

23 

.23 

-n 

-n 

-P 

6 

5.45 

479 

4.79 

+P 

+P 

+n 

0 

0 

0 

0 

-P 

-^P 

+n 

0 

0 

0 

0 

+n 

+  P 

+n 

0 

0 

0 

0 

-n 

+P 

+n 

0 

0 

0 

0 

+P 

-P 

+n 

0 

0 

106 

1.06 

-P 

-P 

+n 

0 

0 

36 

.36 

+n 

-P 

+n 

0 

0 

3 

.03 

-n 

-P 

+n 

0 

0 

0 

0 

+P 

+n 

+n 

1 

.91 

12 

.12 

-P 

+n 

+n 

0 

0 

0 

0 

+n 

+n 

+n 

0 

0 

0 

0 

-n 

+n 

+n 

0 

0 

0 

0 

+P 

-n 

+n 

22 

20.00 

1711 

17.11 

-P 

-n 

+n 

4 

3.63 

584 

5.84 

+n 

~n 

+n 

0 

0 

21 

.21 

-n 

-n 

+n 

0 

0 

5 

.05 

+P 

+  P 
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0 

0 

-P 

+  P 

-n 

0 

0 

0 

0 

+n 

+P 

-n 

0 

0 

19 

.19 

Table    10    (Continued) 
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RAD      NOR      EAS      REC      PREC      SIM      PSIM 


-n 

+P 

-n 

0 

0 

0 

0 

+P 

-P 

-n 

0 

0 

13 

.13 

-P 

-P 

-n 

0 

0 

14 

.14 

+n 

-P 

-n 

19 

17.27 

1617 

16.17 

-n 

-P 

-n 

8 

7.27 

579 

5.79 

+P 

+n 

-n 

0 

0 

0 

0 

-P 

+n 

-n 

0 

0 

0 

0 

+n 

+n 

-n 

0 

0 

0 

0 

-n 

+n 

-n 

0 

0 

0 
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+P 

-n 

-n 

0 

0 

1 

.01 

-P 

-n 

-n 

0 

0 

0 

0 

+n 

-n 

-n 

1 

.91 

97 

.'97 

-n 

-n 

-n 

1 

.91 

69 

.69 

Table  10  (Continued) 
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IV.   DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

The  occurrence  of  earthquakes,  at  the  larger  sample 
sizes,  was  rejected  as  a  uniform  distribution  with  respect 
to  the  lunar  periods,  in  all  cases  except  one  (Table  8). 
On  certain  days  within  the  time  frames  of  SFBAD-I  and 
SFBAD-II,  a  larger  number  of  earthquakes  occurred  relative 
to  other  days  in  these  time  periods.   This  "swarming" 
effect  of  numerous  earthquakes  occurring  within  hours  and 
minutes  of  each  other  may  indicate  the  lack  of  independence 
between  the  events.   The  "swarming"  effect  could  also  indi- 
cate that  earthquakes  do  occur  in  a  non-uniform  manner.   In 
order  to  decrease  the  dependence  between  the  events,  if 
the  dependence  does  exist,  the  data  were  analyzed  according 
to  the  magnitudes  of  the  earthquakes.   This  method  elim- 
inated the  "swarming"  effect.   An  example  of  the  "swarming" 
effect  was  on  June  10,  1969,  when  15  earthquakes  occurred. 
When  magnitudes  greater  than  or  equal  to  2.0  were  examined, 
the  number  of  earthquakes  which  occurred  on  June  10,  1969, 
was  reduced  to  five.   The  null  hypothesis  was  still  rejected 
for  the  three  lunar  periods.   The  fact  remains  that  the 
null  hypothesis  was  rejected  in  all  cases  but  one  at  the 
larger  sample  sizes. 

It  is  interesting  to  compare  the  polar  plots  of  the 
Synodic  Lunar  Periods  of  the  SFVD  and  the  SFBAD-II.   The 
largest  percentage  of  earthquakes  occurred,  in  both  cases, 
in  the  first  interval,  or  at  the  time  of  the  new  moon. 
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Also,  the  time  periods  of  the  two  sets  of  data  were  similar. 
However,  the  polar  plot  of  the  SFBAD-I  was  not  similar,  but 
neither  was  the  time  period. 

The  occurrence  of  earthquakes  with  respect  to  the 
Draconic  Lunar  Period  in  the  SFBAD-I  and  SFBAD-II  are  not 
the  same.   When  analyzing  the  SFBAD-I,  the  null  hypothesis 
was  rejected  at  all  magnitudes,  but  in  analyzing  the 
SFBAD-II  the  null  hypothesis  was  accepted  at  all  magnitudes. 

The  occurrence  of  earthquakes  with  respect  to  the  Anoma- 
listic Month  in  the  SFBAD-I  and  SFBAD-II  are  similar  to 
each  other.   The  rejection  of  the  null  hypothesis  at  the 
one  per  cent  significance  level  was  the  same  for  the  larger 
sample  sizes. 

The  null  hypothesis  was  rejected  v\fhen  analyzing  the 
SFVD  with  respect  to  the  Synodic  Lunar  Period  at  the  one 
per  cent  significance  level  when  using  the  Chi-Square  Test 
but  was  accepted  using  the  Vn  Test  (Both  tests  would  have 
rejected  the  null  hypothesis  at  the  2.5  per  cent  signifi- 
cance level).   The  rejection  was  due  mainly  to  the  occur- 
rence of  17  of  the  110  earthquakes  near  the  time  of  the 
new  moon.   In  comparing  these  17  earthquakes  with  respect 
to  the  combinations  of  the  three  components  of  the  tidal 
force  (as  discussed  in  Section  III),  nothing  of  interest 
seems  to  be  revealed.   Four  of  the  17  earthquakes  had  the 
same  combination  but  the  four  also  occurred  on  the  same 
day.   The  remaining  13  earthquakes  had  different  combinations 
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The  acceptance  in  Section  III,  that  the  distributions 
of  the  110  actual  earthquakes  and  the  10,000  simulated 
earthquakes  were  the  same  seemed  quite  significant.   Further 
analysis  into  why  earthquakes  never  occur  with  certain 
combinations  and  why  large  numbers  of  earthquakes  occur  with 
other  combinations  may  prove  or  disprove  that  the  tidal 
force  can  trigger  earthquakes. 
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V.   THE  COMPUTER  PROGRAM 

The  time  (GMT)  of  each  origin  of  the  lunar  periods 
was  read  into  the  computer  by  the  year  (NA) ,  the  month  (NB) , 
the  day  (NC) ,  the  hour  (ND) ,  the  minute  (NE) ,  and  the 
second  (SE) .   As  discussed  in  Section  II,  the  time  of  each 
origin  was  converted  [Ref.  14]  to  the  corresponding  Julian 
Date  (ORGP(I,J)). 

The  time  (GMT)  of  each  earthquake  was  read  into  the 
computer  by  the  year  (MA) ,  the  month  (MB) ,  the  day  (MC) , 
the  hour  (MD) ,  the  minute  (ME) ,  and  the  second  (SE) .   As 
discussed  in  Section  II,  the  time  of  each  earthquake  was 
converted  [Ref.  14]  to  the  corresponding  Julian  Date 
(JDGMT(I)).   Subsequently  each  of  the  JDGMT(I)  was  con- 
verted to  the  residual  (RESD(I,J))  for  each  lunar  period. 

The  following  list  defines  the  terms  used  in  the  Lunar 
Period  Computer  Program: 

NBOX         The  number  of  equal  intervals  on  the 
interval  (0,1)^ 

0RGP(1,I)     The  origin  of  the  I^^  Synodic  Month. 

0RGP(2,I)     The  origin  of  the  I    Draconic  Month. 

0RGP(3,I)     The  origin  of  the  I    Anomalistic  Month. 

PER(1,I)      The  period  of  the  I^^^  Synodic  Month. 

PER(2,I)      The  period  of  the  I^^  Draconic  Month. 

PER(3,I)      The  period  of  the  I^^  Anomalistic  Month. 

JDAY         Subroutine  to  convert  year,  month  and  day 
to  the  Julian  Day  Number  [Ref.  14]. 
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BEG(l) 

BEG(2) 

BEG  (3) 

EN(1) 

EN(2) 

EN(3) 

JDGMT 
JDLP 

RES 

NCNT 

NPTS 
EXVAL 

PCT 

CHISQ 

SF^ELL 


The  first  origin  of  the  Synodic  Month 

to  occur  within  the  time  frame  of  the  data. 

The  first  origin  of  the  Draconic  Month  to 

occur  within  the  time  frame  of  the  data. 

The  first  origin  of  the  Anomalistic  Month 

to  occur  within  the  time  frame  of  the  data. 

The  last  origin  of  the  Synodic  Month  to 

occur  within  the  time  frame  of  the  data. 

The  last  origin  of  the  Draconic  Month  to 

occur  within  the  time  frame  of  the  data. 

The  last  origin  of  the  Anomalistic  Month 

to  occur  within  the  time  frame  of  the  data. 

The  Julian  Date  each  earthquake  occurred. 

The  time  which  each  earthquake  occurred 

within  the  period. 

The  residual  computed  for  each  earthquake 

and  placed  on  the  unit  circle. 

The  number  of  occurrences  within  each 

subinterval . 

The  total  number  of  earthquakes. 

The  expected  number  of  earthquakes  within 

each  subinterval. 

The  percentage  of  the  number  of  occurrences 

within  each  subinterval. 

The  Chi-Square  statistic. 

Subroutine  to  sort  the  residuals  in 

ascending  order  [Ref .  2] . 
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DIFF         The  Vn  statistic. 

ROOTS        The  Vn  statistic  times  /N. 


47 


0001 


70  00^ 

^  X^  y  sc  V  T- 

c 
c 
c 

C*,*^  tJ-   u.  <Ju  o^ 
T*  *!      -r»  "S^  'f^ 


REAL  ^^8    FRACT,  JDGMT,  JDLPTRES»PER,PERPtnR5,0RGP 

DIMENSION    PER(5  )  t3RG(5  )  ,JDLP(3)  tNCNTOOt  f3j  ,RESD(160,5) 

DIMENSION    CHISQ(5)  ,  SUM ( 5 ) , EM  I N ( 5 )  , EMAX ( 5  J , D I FF ( 5 ) 

DIMENSION    RCOTS( 5) ,0RD(161) 

DIMENSION    ORG P( 3, 10) ,PERP(3,9),PCT(30t5)TEXVAL(5) 

DIMENSION    NUM( 5 )  , JDG MT { 160 )  , BEG ( 3 )  , EN { 3 ) »  ENPTS(5) 

DATA    CHlS0/5*0.0/,  SUM/5^-0.0/ 

DATA    NUM/b-0/ 

RcAD(5,0001)     NBOX 

FORMAT (14) 

BOX^NBOX 

WRITE(6 ,7000)     NBOX 

FORMAT (' 1' ,/,42X, 'UNIFORM    ?• 

<»y  v  V'  Lr  J-  u«  A,  ,v  v  A-  -^  '^-  -v  ^'  4'  -v  •-''  "V  vu  o.'  U'  oj^  %.v  vL*  *.'  «v  "^  ^^ 

<i|^    *^   *>*     p-    *p.    *■•    ^    ij*    A-«    »(•»    Y*    ""™   "T*    ^    TT*    T*   "V    TP    "T*  *"  *  ''*  *;'    1^  ^T"  T*  'l^  *T"  O^ 


ti5 


^    BOXES' J /| 

rf-a    TJS    7C  3|C  3^<k  mpt    tfn  *ft    «^  jJC    ^1^  J^    ?|!    «^   a"^    t 


r    V'     •*—■     ^A*    *^    V*    •** 

■  **■  If*-  ^  T*  nr*  T* 


CONy/ERT    ORIGIN    GREGORIAN    DATES    TO    JULIAN    DATES 

'  T^  1*  '•^  T*  'i*  "P  *t*  "^  *!   "Y*  T*  •**  "C  ^^  ^  ^  ^^  T*  '■^  "*  I*  'i*  V  T^  .    -   n       "  T^  ^'^  ^   •   f^     f^  ^^   ^  ^    '       ^   '^  ^^   'T^  '^^  ^^  ^T*  ^T*  1^  1*  O*  'T' 


6999 


6JJ5 

6000 


6995 
6015 


0040 

0030 


J045 
0035 

0035 

1500 

C»i,  .^  JU  J,  _u 
n^  -r  -nr*  -/-■  V 

C 
C 
C 


0020 


DO    6 

DO  6 
READ 
FORM 
CALL 
ORGP 
ORGP 
CONT 
CONT 
BEG( 
BEG( 
BEG( 
EN(1 
EN(2 
EN  (3 
ORG( 
ORG{ 
DO  6 
READ 
FORM 
CONT 
PER( 
PER{ 
NPTS 
DO  0 
DO  0 
NCNT 
CONT 
CCNT 
DO  0 
DO  0 
RESD 
CONT 
CONT 
11=1 
READ 
FORM 
IF(M 


300 
005 
(5, 
AT{ 
JO 
(I, 

(It 

INU 
INU 

1)  = 

2)  = 

3)  = 
)=C 
)=0 
)=0 

4)  = 

5)  = 
015 
(5, 
AT( 
INU 
4)  = 
5J  = 
=0 
333 
040 
(I  , 
INU 
INJ 
33  5 
045 
(I, 
INU 
INU 


1  =  1,3 
J=l,13 

6999)     MA, 
IX,  14, I  3, 


NB 


y  NC  , 
t  14, 


ND,NE,SE 

I3,F6.2) 
AY(NC,M3,NA, JD) 

J)=JD+(FLDAT{ND)+FLOAT(NE) /60. 0+SE/ 3600.0 ) /24.0 
J )=ORGP(  I , J  )-0.50 


0RGP{1,2)     . 
0RGP(2,2) 
0RGP( 3, 2) 
RGP(1,3  ) 
RGP(2,9) 
RGP{3,9 ) 
0RGP{2, 1) 
OR&P( 2, 1) 

1  =  1,3 
6995)     { PERP(  I , J)  ,  J  =  l,9) 
10F3.5) 

27. 32156 
27.32158 

1=1, NBOX 

J=l,  5 
J)=0 
E 
E 

1=1, lb3 

J=l,  5 
J)=0.0 
E 
E 


( 5 , 1 5  03 , EN  D  =  302  3 )    MA , MB , MC , MD , ME , S  E 
AT(  IX,  14,1  3,  13,  14,  13, "^6. 2) 

A. EQ.O )    GO    TO    0020 

•"    1"   *•    T*  "^  1*  -■-    *r*      *    V*     *    -*    *r-   -p    -t*.   —    -;■»   -f^    -f.   ^^    TT-    ^    .,-    T^    ^p   <!-.    -,^    ^    -*    ^-    ,^   -|t  ^   .y.    -X-    *?    *r*  T  • 


•     ^    *^    OL     <l-     V^     ■A'    U^     «,/«    W'    oA^ 


CONVERT     EARTHQUAKE    GREGORIAN    DATES    TO    JULIAM    DATES 


*^  Of  •^  U|^  >t^   J* 

CALL  J 
NPTS=N 
FRACT= 
JDGMT{ 
GO  TO 
DO  430 
DC  400 
IF( I.G 
I  F  (  J  DG 
IF( JDG 
DO    603 


^-  ^  ^^  ^ 

DAY( 
PTS  + 
(  FLG 
NPTS 
0005 
0  1  = 
5  L  = 
T.3) 
MT(L 
MT(L 
0    J  = 


MC,MB,MA, JD) 

1 

AT( ^D ) +F L 0 AT ( ME )/60. 0+SE/ 3600.0) /24.0 

)=J0+FRaCT-0.500 

1,5 

1,NPTS 
GO     TO    4100 

)  .LT  .3EG{  I  )  )    GO    TO    4005 
) .GT.EN(  I  ) )     GO    TO    4335 
2,10 
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IF( JDGMT (L ) .Cr  .DRGP(  I  ,  J  J  )     GO    TO    6030 

ORGd  )  =ORGP(  I,  J-1  ) 

PER( I )=PERP(I, J-1) 

J=10 
6030    CONTINUE 
4100    NUM(  I)  =  MUM( I  )  +  l 

JOLP{ I )=JDGMT(L)-CRG(I ) 

f      «y«  *j*  *  •  3|*  9p  ^p-  *(»  ^f*  *^  JsC  ^«  Jp  *  ■«  rf^  X|C  SJC  ^C  jrT    3jC  ^  9C  SJC  SC  ^jp  «f«  3fC  #^  Jya  ^»  ^»  jJC  3|t  "^  ^i-k  *j^  ^ .  *r»  *{*  *if  3^^  3p  ^f*  •  *  i^  •!*  ^"  *t*>  *y~  n^   "t*  »y*  •*'^  'i*  *■  -  ^t"  "  ■*  "^  t*  X* 

c 

C  COMPUTE    RESIDUALS 

C 

RES  =  Df10D(  JDLP{  I),PER(I))/PER(I) 
IF(RES.  LT.  3.CO)     RES=RES+1.D0 

C^U  *Xr    O^  «■#   «  «     aJU  ^^    ^^    '■'^     "V    ••><'    0-<     «U    ■hL'   ^<-    V.U    aXf    ij.      Jl.r     *.*•    Ov    «^    >iJ^    ^<    .jU  O-r    ^y    iJU    •*l'     a^    vV    Jl'    «<'    Ov    -^    •''    aJU   (.U     •A'     fciU    w^    «!•    JU     M     J^    iJU    OU    iJU    ^V    *X.    <.t«    vlv    J<     *^     •*-    ^     •fU     JU    JU 
-1*  '.*  V  '.*  '  '   ^  -T^   "T^  *.'  'T'  '.*  ■"   T*  "V-  n*  ':-   -^  •^-  ¥^    T-  -  -  -^  'r  "^  'T  -1'  *]*  -^  *«"   -r-  **  'f-   '*   '-^  ',"  -/*  T"   -.^   1'   -T  ^   ^   -T   1* 

c 

C  PLACE    RESIDUALS    INTO    20    BOXES    ON    INTERVAL     (0,1) 

C  FOR    EACH    LUNAR    PERIOD 

C 

INDEX=DES--;=BGXf  1.0 

NCNTdNDEXt  I)=NCNT  (  INDEX,  I  )+l 

RtSD( I  I , I )=RES 

II=I I+l 
4005    CONTINUE 

ENPTS( I )=NUM( I ) 

EXVAL(  I  )  =  ENPTS ( I )/BGX 

DO    0070    J=1,N3  0X 

PCT{  Jt  I  )=FLOAr  (NCNTC  J,  U  )  /ENPTSd) 
0070    CONTINUE 

C 

C  COMPUTE    CHI-SQUARE    ON    SOXES 

C 

CV'*'^'  «^  %U  •.  .•  k>«  «.  .  J^  *i«  ijf-  OU  0<  ^r  -vV  o^  >u  <>;^  *,V  V'  •*'  ^''  •^'  •V  '•'^  -kJ^  "i^  <^'  ■^''  V' 
^*.tf^  -^o  T-  '*'  •^'  -i*  -r*  "  *■,*.  iA>  ^  .  -7^  -,-  .^«.  *i(~  -,-  .-y^  ^.  *,.  jr*  ,^  .:^  .ft  rf  »i-  *■,*  «■  n  *i-  f-   '•'*  •^"'  i^  •'r*  ^  '\'  T-  T  *i*  -r-  ^-^  t^  *■'  ^  */*  *-p  -y.  ^*  -y.  ^  -,»  Jj»  -i'.  *v  '■.■'  t-  ».*  t* 

DO    0098     J=ltMBOX 

SUM(  I)  =  (  FLOAT  (  NCNT(  J  ,1  )  )-EXVAL(I)  J  =?-2/ EXV  AL  (  I  ) 
CHISQ(  I  )=CHIS3(  I  )+SUM(  I ) 
0098    CONTINUE 

*  ■>■  *Y~  ""I*  'T^  ',~  T*  'i"  'V^    ^r*  'f   "'•  T*   'P  n*  *"*  T"  "<^  1^  "^  ^■^  IT*  ^  If  "^^  n^  n   T^  T'  'V^  T*  n'*  nf*'  ^^  T"  1*  f*'  t"  T^  T*  i"   ^*  *1 '  nr*  "p  1^  1^  "  "T*  "T"  ^r"  ^"^   **  *"   •*r^  "IP  I"  ',•  *!* 

c 

C  SORT    RESIDUALS 

C 

CALL    SHELL(NUM{  I  ),  RESDd,  n  ) 
NO=NUM( I ) 
DO    0033    J=1,N0 
ORD{ J)=FL0AT{J-1)/ENPTS(I) 
0033    CONTINUE 

ORD(NJM{ I)+l)=1.0 
EMINd  )  =0.  0 
EMAX( I )=0,0 

C 

C  COMPUTE    MODIFIED    K-S    TEST    ON    THE    LUNAR    PERIODS 

C 

■n^  T"  ',"'  V  <-v  *p  -™  -|%  ^  ',*  -n*  V-  -T«  --  -T*  ■»[*■  -^  jx  ^  -p  -.*  ^~  ^  .-I*  .  .  *t.  * ,  .■r*  ^  «-«  '^  '^r'  'r  "^  ^  'T'  '<''  ^  T-  -v  'r  t  -t*  i*  nr  -  -  -tt  -t*  t  t*  t*  '^  -*  ^  ■'^  '  ^  'k-  n~  "* 

DO    0054    J=1,N0 

RESDX=RESD(J, I ) 

EMAX( I } =AMAX1( RESDX-ORO( J) tEMAXCI ) ) 

EMINd  )  =  AM  INK  RE  SO  X-ORD(  J  + 1  )  ,  EMIN(  I  )  ) 
00  54    CONTINUE 

DIFF{ I )=EMAX( I )-EMIN( I ) 

ROOTS (  I  )=DIFF(  I)vSQRT( ENPTS(  1}  ) 

11=1 
40  00    CONTINUE 

«.  *,.  ,v*  .,,  r:  ^  „  J^i  X-.   j^  ,r-  1*  1^  -■  ■  ^  '.-  '-  -r^  *F  »,*  -^  -■  -  -Y^  -r*  •■ '  -v-  ^^  '.*  T-  -f»  1-  -!■•  '^'  ^•.  **  3p  3X  ^Y*  '.•*  -t~  •■•■  r'^  'r'  -r-  -r>  -r-  -^  -^  -r-  or-  -^  **  -r  "T*  »•."  -r  n*  'T*  ^ 

C 

C  TABULATE    THE    RESULTS 

C 

WRITE! 6,9999) 
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9999     FORMAT (9X, • SYM3DIC    MON TH ' , 5X, ' DRACON I C    MQNTH«,5X, 
1' ANOMALISTIC    MONTH • , 5 X , ' S I DER EAL    M0NTH',5Xt 
2'TROPICAL    MONTH' ,//) 
DO    0060     I=1,NB0X 

WRITE (6, 9903 )(  NC NT (  I, J)  ,  J  =  l ,5) 
9900     FORMATC  I5X,I4,li>XT  I4,18X,I4,15X,I4,15X»I4)     ' 
0060    CONTINUE 

WRITE  (  6,4999)  (  NUM(  J  J  ,J=1  ,!5) 
4999    FORMAT (/,11X, I  4,  '     QUAK E S • , 8X, I 4t '     QUAKE S •  ,  1 1 X  ,  I  4 , 
1'     QUAKES' ,8X, 14, •     QUAK ES ' , 8X , 14 , •     QUAKES',/) 
WRITEC  6, 7005) 
70  05    FORMAT (//,4oX, 'PERCENTAGES' ,/) 
WRITE (6,9999) 
DOO155I=l,N30X 
WRITE(6,8100)(PCT{ I,J) ,J=1,5) 
8100    FGRMAT(I1X,F7.3,12X,F7.3,12X,F7.3,15X,F7.3, 12X,F7.3) 
0155    CONTINUE 

WRITE(b,2044) 
2044    FORMAK • 1' ,// ,45X, 'CHI-SQUARE    TEST',/) 
WRITE(6,9999) 

WRITE(6 ,32  00)( CHISQ(J) ,J=1,5) 
82  00  FORMAT(  11X,F7. 3,12X,F7.3,12X,F7.3,15X,F7.3  ,12X, 
1F7. 3, //////////) 
WRITE(6,9999) 

WRITE ( 6,4988)  { NUM{ J )  ,J=1,5) 
4988  F0RMAT(2X,9X, 14,  '  QUAK ES ' , 8X , I  4,  '  QUAKE S '  , IIX , I  4 , 
1'  QUAKES' ,8X, 14, '  QUAK ES ' , 6X , 14 , '  QUAKES',//) 
WRITE (6, 2022)1 DI^^CJ)  ,J  =  1,5) 
2022  FORMAT (IX, 'MODIFIED',/ ,3X,  'K-S •,3X,F 9.5, 10X,F 9. 5,1  OX, 
iF9.5,13X,F9.5,13X,F9.5,/,lX,'STATISTIC',//) 
WRITE(6,2033)( ROOTS (J ) ,J=1,5) 
20  33  FORMAT (2X, 'ABOVE', /,1X,'STATISTIC',/,2X,«TIMES',2X, 
lF9.5,10X,F9.5,10X,F9.5,i3X,F9.5,10X,F9.5,/,lX, 
2  'ROOT(r;PTS)  '  ) 
STOP 
END 

SUBROUTINE  jrAY( D,M,Y, J ) 

INTEGER ^4  D,Y,YA,C 

IF(M.GT.2)  GO  TO  0055 

IM=M+9 

IY=Y-1 

GO  TO  0065 
0055  IM=M-3 

IY  =  Y 
0065  C=IY/100 

YA=IY-100"C 

J=( 1460  97-C) /4+( 1461^YA) /4+ ( 153- I M+2 ) /5 +D-H72 1 1 1 9 

RETURN 

END 

SUBROUTINE    SHELL(N,X) 

DIMENSION    X(l) 

M  =  tsl 
0101    M=M/2 

IF{M.EQ.O)     RETURN 

K  =  N-M 

J=l 
0111  I=J 
0121  IM=I+M 

IF(X(  I  )  .LE.Xd  M)  )  GO  TO  0131 

S=X(I) 

X(I) =X( IM) 

X  (  I  M  )  =  S 

I=I-M 

IF( I.GE.l)  GO  TO  0121 
0131  J=J+1 

IF(J.LE.K)  CO  TO  0111 

GO  TO  0101 

END 
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